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Abstract 

By presenting a relation between average energy of the ensemble of probe photons and energy 
density of the Universe, in the context of gravity’s rainbow or doubly general relativity scenario, 
we introduce a rainbow FRW Universe model. By analyzing the fixed points in flat FRW model 
modified by two well known rainbow functions, we find that the finite time singularity avoidance 
(i.e. Big-Bang) may still remain as a problem. Then, we follow the “Emergent Universe” scenario 
in which there is no beginning of time and consequently there is no Big-Bang singularity. Moreover, 
we study the impact of a high energy quantum gravity modifications related to the gravity’s rainbow 
on the stability conditions of an “Einstein static Universe” (ESU). We find that independent of a 
particular rainbow function, the positive energy condition dictates a positive spatial curvature for the 
Universe. In fact, without raising a nonphysical energy condition in the quantum gravity regimes, 
we can address an agreement between gravity’s rainbow scenario and basic assumption of modern 
version of “Emergent Universe”. We show that in the absence and presence of an energy-dependent 
cosmological constant A(e), a stable Einstein static solution is available versus the homogeneous and 
linear scalar perturbations under the variety of obtained conditions. Also, we explore the stability 
of ESU against the vector and tensor perturbations. 
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1 Introduction 


In the framework of general relativity (GR), the gravitational force is explained in terms of the space- 
time curvature so that the field equations connect the space-time geometry to the matter content. More 
technically, GR exhibits a Universe modeled by space-time with a mathematical structure formed by four 
dimensional differentiable manifold [T] . It is commonly believed that a unique mathematical framework 
to reconcile the quantum mechanics (QM) with GR is highly dependent on our understanding of the 
space-time geometry. In other words, there is a comprehensive agreement that the geometry of space- 
time is fundamentally explained by a quantum theory. It is assumed that the Planck energy Cpi addresses 
a critical energy scale of transition from classical GR to quantum gravity (QG), a theory which is 
supposed to challenge the most fundamental and long standing issues of modern physics. Despite the lack 
of a complete theory of QG, it seems a semi-classical or effective phenomenological approach of QG may 
guide us to disclose the mysterious nature of QG [U [3111 [5]. It is important to understand how to extract 
the testable predictions from this fundamental theory. Interestingly, all semi-classical approaches offered 
so far unanimously insist on the existence of a minimal measurement length in the nature, known as the 
Planck length Ip. Besides, it is believed that the modified mass-shell condition (or dispersion relations) 
arising from deviation of Lorentz invariance in effective phenomenological approaches HIS], may justify 
some of the phenomena taking place in astronomical and cosmological scale, like the threshold anomalies 
of ultra high energy cosmic rays and TeV photons [iiz]. It is noticeable that such Lorentz invariance 
violation is already predicted in the context of other approaches to QG such as: non-commutative 
geometry [7], spin network in Loop quantum gravity [8], and string field theory [9]. Throughout the 
present work, our attention specifically is focused on a semi-classical formalism knows as “Gravity’s 
Rainbow” which has been designed by Magueijo and Smolin m- In fact, it can be said that the 
rainbow gravity is nothing but the doubly special relativity (DSR) [3] in the presence of curved space- 
time which is knows as the “Doubly General Relativity ”(DGK). In this formalism, there is no single fixed 
space-time background, namely the space-time background appears as a geometric spectrum in terms 
of the energy scale of the particle probe. Therefore, in the cosmological setting, we are dealing with the 
rainbow modified metric by introducing a function in the Friedmann-Robertson-Walker (FWR) metric 
which depends on such variable energy scale. Formation of the rainbow metric proposal involve various 
motivations such as the lack of a trivial definition of the dual position space in DSR [10]. However, this 
approach to QG is not problem free and considerable number of works have been focused on its various 
aspects m- Somehow, the proposal of gravity’s rainbow is similar to the idea of “Running Goupling 
Gonstants” in particle physics and field theory, in the sense that at very small distances (high energy) the 
space-time geometry is related to the energy scale of the particle probe m- It is worth mentioning that, 
apart from the effective approaches to QG like DSR and DGR, it has been argued that if the physical 
space-time of the standard GR theory is emergent, then it is expected that we face with a radical picture 
of the Universe at fundamental quantum scale m- The basic idea of “Emergent Universe” dates back 
to the works of Eddington in 1930 m, inspired by the proposal of Einstein static Universe (ESU). 
In simple words, ESU suggests that the content of the Universe as a closed system is made up of at 
least cosmological constant and normal matter. While ESU in the context of GR is not stable versus 
spatially homogeneous perturbations, it predicts that the Universe in past might have been emerged as 
an static initial states [14]. We know that the data obtained by cosmic microwave background (GMB) 
observations, generally supports the inflationary scenario. Also, by implementation of some singularity 
theorems with the geometrical assumptions that AT < 0 or iF = 0, it has been shown that despite the 
occurrence of inflation, the Universe had a violent beginning in the past m- Hence, one of the most 
important motivations for studying QG is to avoid the singularity at the beginning of the Universe. 
Recovery of ESU in the framework of inflationary Universe by Ellise and Maartens, knows as modern 
version of Eddington emergent Universe, was in this direction of research interest m- According to 
the modern standard cosmology, the inflationary expansion of the Universe has quickly eliminated any 
original spatial curvature. On the other hand, the GMB observations implies Uq > 1 (1.02±0.02) for the 
total density parameter. This means that the geometry of Universe is not exactly flat, rather it could 
have non-zero spatial curvature at the beginning, resulting in negligible late-time effects. In this sense, 


2 


the emergent Universe scenario allowed the existence of a positive curvature Universe that emanates 
asymptotically as a static initial state known as ESU which afterwards experienced the inflation. The 
stable ESU addresses the existence of a fixed point around which our Universe in the early times was 
eternally fluctuating. In the case of the Universe filled with a massless scalar field (p with a suitable 
inflaton potential, these fluctuations had been disappeared and the stable state of Universe had been 
turned to inflationary phase. Within the context of standard cosmology, the prerequisites of a one 
or further flat wings (of course with a little positive gradient) in the inflaton potential, can make an 
emergent Universe. This transmission from a stable ESU to inflationary phase takes place around the 
different values of e-folds. For instance, in m by analyzing the spectrum of inflationary perturbations 
of some stable cosmological models in the context of GR, it has been demonstrated that the transmission 
occurs in over 60 e-folds. Also, it should be noted that while these inflationary emergence GR-based 
models respect to CMB constraints, they suffer from a fine-tuning problem. This problem can be 
resolved in the context of modified emergence models of GR. However, recently in m some arguments 
are presented which tries to show the emergent cosmological scenario can not really be past-eternal. As 
a prominent feature of emergent Universe scenario, one can point to the removal of initial singularity 
and also the horizon problem. This is a strong motivation to study the ESU models along with their 
stability conditions in the presence of high energy corrections of GR. For instance, one can point to 
the works done in the context of massive gravity |18j , Hof ava-Lifshitz model of gravity m, braneworld 
scenarios [50], induced matter theory loop quantum cosmology [55], /(A), f{T) and /(G) gravity 
[50] . As discussed at the beginning of this section, gravity’s rainbow or DGR is also considered as 
an alternative of GR with high energy corrections. Besides, it has been shown that in the study of 
FRW cosmology and isotropic quantum cosmological perfect fluid model in rainbow gravity setup, some 
conditions are derived which prevent the initial singularity [25] . Inspired by this introduction, our main 
goal in this paper is to study the stability of ESU against the homogeneous scalar, vector and tensor 
perturbations in the presence of rainbow’s gravity setup as a QG modification of GR. But before doing 
this, by following the phase space method we will analyze the status of the Big-Bang singularity in the 
framework of a flat rainbow FRW Universe model. 

The present paper is arranged as follows: In section 2, we have derive rainbow modified Friedman 
equations in details for a general non-flat universe. In section 3, by introducing two different forms of 
rainbow functions depending on energy density, firstly in absence of cosmological constant, we obtain 
the stability conditions of ESU against the homogeneous linear scalar perturbations. In the follow of 
this section, we pursue our goal in the presence of a cosmological constant related to energy density via 
introducing a new rainbow function. In section 4, we analysis the stability against the vector and tensor 
perturbations. Finally in section 5, we give the summery and conclusions. 


2 Gravity’s Rainbow Modified FRW Universe 


In this section, firstly we have a quick review on the gravity’s rainbow theory, known also as DGR. 
Then, with note to the importance of spatial curvature in Emergent Universe scenario, we derive the 
modified Friedman equations for general non-flat universes. In DSR, the modified dispersion relations 
of a massive particle with mass m reads as 

eV^(e/ep,??) -/ff^(e/ep,?7) = , (1) 


where f{e/ep,r]) and g{e/ep,r/) are the well known energy dependent rainbow functions and 77 is a 
dimensionless parameter. In the low energy limit, the modified dispersion relation (HD reduces to the 
relativistic dispersion relations. So, the rainbow functions f{e/ep,r]) and g{e/ep,r/) satisfy the following 
conditions 


lim f{e/ep,r]) = 1, lim ^(e/ep, 77 ) = 1. 

e/cp-^O e/cp-^O 


( 2 ) 


Indeed, the condition ([5]) points to the correspondence principle. Based on the discussion carried out 
in [^, for established position space in DSR, theories including free fields must also lead to the plane 
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wave solutions in flat space-time, despite satisfying the modified dispersion relations ©• For this reason, 
contraction between infinitesimal displacement dx°‘ and momentum pa, must be linearly invariant, i.e 

dx°'pa = dte + dx^pi . (3) 


In fact, the linear contraction Q guarantees the existence of the plane wave solutions. The rainbow 
metric generally has the form of 


ds"^ 


-1 


/^(e/ep:?7) 


dt^ 




dx^ 


(4) 


Using a one parameter family of energy momentum tensors, the gravity’s rainbow modified Einstein 
equations will be written as 


Ga/3(e/ep) = S'KG{e/ep)Tap{e/ep) -b ga/3(e/ep)A(e/ep) , (5) 

so that G{e/ep) = hi{e/ep)G and A{e/ep) = /i 2 (e/ep)A where/ii(e/ep) and /i 2 (e/ep) are energy dependent 
rainbow functions. This means that in DGR setup, the Newtonian gravitational constant and the 
cosmological constants are energy dependent such that in low energy limit, we recover G(e/ep) = G and 
A(e/ep) = A. In order to achieve our main goal, in this section we will set A = 0 for simplicity, and 
begin with the following modified FRW metric of a homogeneous and isotropic Universe 
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where hij represents the spatial part of the metric. Commonly, it is assumed that the energy e of particle 
probe for each measurement is constant and independent of space-time coordinates. Nevertheless, such 
assumption for the measurements at early Universe seems to be far from reality. Therefore, it is expected 
that the background metric of space-time throughout its evolution is affected by the energy e of particle 
probe. Then, it is reasonable to consider the evolution of energy e with the cosmological time, denoted 
as e{t). In this regard, we will derive the modified FRW equations. In doing so, the following ansatz is 
usually suggested m 

= (7) 

where / has the form of f{e/ep). Therefore, the modified Einstein field equation can be written as 

Rai3 = -8TTG{e)Sap{e), ( 8 ) 


where Rajs is Ricci tensor dehned as follows 


= ^ - ^ + r^^F^, - F^^F^^, a, /3 = 0,1, 2, 3, 
and Sapie) is written in terms of the energy momentum tensor Tap{e) as 

Sapip) = Tap{e) — -^gapTp{€) . 

To continue our calculations, we need the non-zero components of the affine connection as 
r!So = -p FO. = Aaa%, z,j = 1,2,3. 


(9) 


( 10 ) 


( 11 ) 


It is seen that unlike the usual FRW metric, for modified FRW metric ([6]), the components of the affine 
connection with two time indices remain non-zero. By putting m into and after a straightforward 
calculation, we obtain 


Rqo — 3— + 3—^ , 
a a} 


( 12 ) 
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and 


(13) 


Rij = Rij - ((ad + 2b?)- ffaaj 6ij , 
where Rij denotes the purely spatial Ricci tensor defined as 

D. _ m _ I -pm-pm _-pmpn 

~ dxi dx^ ^ ^ ^ ' 


(14) 


The spatial components F® j, of the four dimensional affine connection are identical with those of affine 
connection computed in three dimensions from the spatial three-metric i.e. 



/ dhjji 

\ dx’^ 


dx^ 


^^jk \ _ pi 

dx^ ) ~ ’ 


(15) 


where F* denotes purely spatial affine connection and is the inverse of the 3x3 matrix hij. Therefore, 
we have F^- = kx^hij which results in 

Rij = —2khij . (16) 

Here, k has a geometrical interpretation. Indeed, it measures the spatial curvature with zero, negative 
and positive k values corresponding to flat, open and closed Universes respectively [5S]. The Ricci tensor 
expression m can be rewritten as 


Rij = — ((ad -I- 2a^)f^ — ffaa + 2k^ hij . (17) 

Now, for obtaining the components of SajS, we consider the perfect fluid with the energy-momentum 
tensor as 

Tap = (p + p)UaUp + pgaP , (18) 

so that p and p represent the energy density and the pressure, respectively. Also, Ua is the velocity four 
vector defined by Ua = (/“^, 0, 0, 0), with the norm of g°‘^UaUp = —1. We should note that the diagonal 
components of the modified FRW metric tensor gap is {—f~‘^,hii) with the signature (—, -I-,-I-,-I-) so 
that the components of spatial tensor ha are independent of the rainbow function /. Using the equation 
one can decompose the time and spatial components of Sap tensor as follows 



Stt — Ttt -|- —gttT, Sij — Tij — —hija?T , 

(19) 

where 

Ttt = pf~'^, Tij = hijPp , 

(20) 

and 

T = g^^Tap = {-pTip). 

(21) 

By substituting expressions (1^ and (1^ into (ITOl). we have 



Stt = ^{p + ip), S^j = ]^{p-p)Phij. 

(22) 


Ultimately, the first and second rainbow modified Friedman equations for the metric parameterized by 
the varying energy probe ([6]), takes the form of 


^a^ a2/2 3 /2 ’ 

d 47rGp(l -I- 3aj) 1 a f 

a 3 P a f ’ 


(23) 

(24) 
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Here, also for simplicity, it is assumed that gravitational constant G is independent of energy e . It 

should be noted that the existence of / term in the second rainbow modified Friedman equations, does 

not mean the explicit dependence of rainbow function to the time i.e. /(t), at all. As explained above, 

/ is an explicit function of the energy of the test particles which are probing the geometry of space-time 

in early Universe. Given the fact that e can vary with respect to the evolution of Universe, so / can be 

an implicit function of cosmic time and not explicit. We note that, the explicit dependence of / to the 

time, may be leads to this wrong result from equations (I23|) and (l24l) that they are nothing but the usual 

Friedman equations of GR so that the rainbow function / in this case is merely the gauge parameter 

determining the choice of time. Hence, based on the gauge freedom, one may choose the gauge of / = 1. 

In contrast to this misconception, the gravity’s rainbow theory is a high energy modified theory of GR 

in which according to correspondence principle, for the case of low energy limit i.e. ^ ^ 0 (/ —>■ I), 

modified Friedman equations (1^ and (IMl) reduce to the standard equations for the FRW universe. 

Also, we mention that due to the existence of the rainbow functions, t may play the role of proper time 

based on the gauge choice N = [53]. As expected, in the limit of GR, t represents the proper 

J2 

time indicated by the choice of A^ = Also, by combining the modified Friedman equations (E31) and 
([M|. we get the following energy conservation equation 

p -\- 3—p(l -|- u)) = 0 . (25) 

a 

It can be mentioned that in framework of DGR, the form of equation of state (EOS) p = ujp, remains 
unchanged for massless prob particles such as photons while for the massive ones, this equation of state 
will be modified, see [29l [10] for more discussion. 

3 Einstein Static Universe and Scalar Perturbations 
3.1 In the Absence of Cosmological Constant A 

In this section, we plan to apply the linear homogeneous scalar perturbations in the vicinity of the 
Einstein static Universe and explore its stability against these perturbations. To achieve our goal, we 
need to fix the rainbow function /. To this end, using dispersion relation offered in [30] . we introduce 
rainbow functions as follows 

/(e) = (l-—)-i, (26) 

epz 

where via suggesting the average energy e = (c is some constant) [31) . this rainbow function will 

be rewritten as 

f{p) = {l-\ip-^)-\ (27) 

where ^ The authors of [3T], to reach the above relation considered a large ensemble of prob 

photons which are in thermal equilibrium. This assumption is reasonable since based on the standard 
model of cosmology, early universe passed a radiation dominated era with p = ^p. One may ask the 
question that why the identification e ~ e is used to obtain the rainbow function (EZI) ?. The answer is 
that, we are deal with the energy of prob particle in the rainbow metric as the statistical mean value of 
all prob photons in radiation domination. Indeed, we deal with the average effect of photon particles in 
radiation dominated era and no with a special elected photon from the radiation m- It is worthwhile 
to remind that the form of e in terms of p is independent of the modified dispersion relation picked for a 
specific model [23]. We also mention that the above rainbow function is valuable from the theoretical 
viewpoint so that in [30] it is shown that in the absence of the varying speed of light (VSL) proposal, it 
can removes the horizon problem in early Universe. Moreover, it is seen that when p or e get their largest 
values, then /(p) or /(e) becomes infinite which causes the time-like component of the rainbow metric 
vanishing. This means that the FRW metric which is modified with rainbow function (1271) becomes 
degenerate at that time and has no inverse. Indeed, a degenerate metric addresses the existence of 
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other distinct possibility for lightlike dimension (other than timelike and spacelike dimensions). Also, 
remember that in the Palatini formulation of standard GR, degenerate metric also appears. One of 
the consequences the degenerate metric is that the curvature remains bounded and the topology of 
space-time can change. Overall, it is believed that singularities arisen from degenerate metric have a 
milder manner than other type of singularities and seem appropriate for a QG proposal [32] . Here, it is 
necessary to review the finite time singularity issue, including the Big-Bang singularity, in the presence 
of rainbow function (1271) . Indeed, we want to investigate whether the presence of rainbow function (1271) 
in this system will lead to resolve the singularity problem. Inspired by the idea proposed in [31 [SI], 
we want to resolve this problem by finding an upper bound on the density of energy p. In other words, 
according to I31I3I] the finite time singularity issue will be eliminated via the presentation of a fixed 
point as pf for energy density p which is reached at an infinite time. More exactly, the author of [31j . by 
following the terminology used for stability analyzing the dynamical systems in [33] , demonstrated that 
for any first order system as p = 0{p), the finite time singularity will be solved through the fulfillment 
of either of the following conditions: 1) the function 0{p) be a continuous and differentiable on a range 
enclosed by zeroes of function 0{p). 2) asymptotically, function 0{p) growths like a linear function as 
K{p) or slower than it, i.e K{p) > 0{p). Therefore, for cancelation of the finite time singularity, it is 
sufficient that one of these conditions is satisfied. We begin our analysis in this way by substituting 
the rainbow function (1271) into energy conservation equation (I25|) and modified first Friedmann equation 
(ESI) (for case fc = 0) to obtain the ordinary differential equation (ODE) as p = 0{p) where 

0(p) = -4p(l-Cp^)(^—pj . (28) 

The advantage of the dynamical system analysis method followed in isaisT] is that by having fixed 
points and regarding the asymptotic behavior of 0{p), one is able to predict the demeanor of the system 
without need to have a detailed form of the solutions. Hence, by setting w = ^ (since our attention 
is on radiation dominated state in order to study the initial singularity), equation ((28l) results in the 
following two fixed points 

Ph =0, Ph = C-ff- (29) 

However, to understand the qualitative manner of a solution, we should know that how long it takes to 
get a fixed point. By a straightforward calculation, one can show that the time required to reach these 
fixed points, can be obtained as 


f=_r— ^ (30) 

where p* represents an arbitrary initial finite value for density which lies in the intervals between the 
fixed points. The negative sign in the back of integral relation, denotes a backward in direction of time. 
We find that for the case of pf^ = 0, the integral (ISOl) does not converge i.e. |t| —)> oo, while for the fixed 
point pfy = fbe integral (ISOl) is not solvable and its numerical values fail to converge. At this 

point, let us to draw the plot 0{p) — p (or p — p) in the figure 1, in order to get a qualitative analysis 
of the situation. As it is seen from the figure 1, 0(p) has two fixed points = 0 and p/j = 1. In 
fact, any of these fixed points are equivalent to a de Sitter space since in these points we are dealing 
with constant solution (p = 0). The first p/j is a future fixed point because for the case p* > 0 we have 
0(p) < 0, while the second pf^ is an early (or past) fixed point since for case of p* < 1 and p* > 1 we 
have 0{p) <0 and 0(p) >0, respectively. It is noteworthy that for distinction between the type of this 
singularities, we follow Ref. [31]. These fixed points classify the possible solutions into two classes: 1) a 
solution belongs to interval p G [0,1], if p/^ < p* < pf^ and 2) a solution belongs to interval p G [1, oo), 
if p* > pf^. Then, to get the fixed point pj^ by beginning from some initial value as p* and also from 
p /2 to p*, the time required is infinite. In this interval, 0(p) is continuous and differentiable. Then, 
according to approach in m by pursuing the terminology of |S3], one can say that the first solution 
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0.2 


Figure 1: The behavior of 0(p) (1281) in terms of p for flat rainbow FRW Universe model. To simplify, we set 
epi = c and 8ttG = 3. Fixed points are located in p/j = 0 and pf^ = 1. 


is free of physical singularity so that it interpolates steadily from p/^ to p/j. It is clear that the first 
condition already is not usable for the second solution, which represents a clear violation of the second 
condition. On the other hand, for the case p* > p/j, the function (l28ll is growing quicker than a linear 
function. So, we conclude that the second solution is not free of a finite time singularity. Because we 
are interested in studying initial singularity problem, so let us mention to stability status of early fixed 
point p/j. A fixed point is stable if by putting a neighboring initial value, the trajectory of the solution 
remains always near to the fixed point. Equivalently, a fixed point is unstable if for any point at vicinity 
of the fixed point, one can find some solution that starts near the fixed point but go away from it in a 
finite time. The fixed point p/^ is an unstable point because Ip=p /2 >00- Then, although it takes 
an infinite time to get a fixed point, the issue of finite time singularities avoidance remains unsolved and 
eventually this fixed point will collapse. 


Now, we investigate the DGR theory form the emergent universe point of view in which there is 
no beginning of time and consequently, there is no big bang singularity in the early Universe. In this 
scenario, the Universe did not born from a Big-Bang singularity in a past finite time and rather it 
possesses an eternal Einstein static state. In this context, the key point is the required conditions for 
the stability of the existing ESU, in which we will explore in the following of the present paper. The 
ESU in the DGR scenario with rainbow functions varying with cosmological time can be obtained by 
the conditions d = d = 0 through the equations (1251) and (125)) as 


SttGpo 


(31) 


and 


47rGpo(l + Sui) 

3/(po)^ 


= 0 


(32) 


where ao, po and oj denote the scale factor, the energy density of the Einstein static Universe and 
barotropic equation of state parameter (p = cap), respectively. Considering the positive energy condition 
Po > 0 through the equation it is seen that for the Einstein static Universe, the spatial curvature of 
the Universe should be positive. A: > 0. Also, Eq.([32|), shows that for the Einstein static Universe in the 


^At a fixed point where 0{pf^) = 0, if \p=pf^ > 0, we have increasing O at pf^, or equivalently /(p /2 — 5) < 0 < 

/(p /2 + (5) for all sufficiently small and positive 5. This means that if we start with initial value p* > p /2 th® vicinity 
of p/ 2 , since /(p*) > 0 then the trajectory of the ODE solution increases its value of p /2 and moves away from the fixed 
point. Also if we start with p* < p/ 2 , but near to p/ 2 , again trajectory of the ODE solution moves away from the fixed 
point but by decreasing its value of pf^. Therefore, if ~^~\p^pf 2 ^ conclude that the fixed point pf^ is unstable. 

Similarly, one can show that if \p=pf^ < 0, the fixed point will be stable. 











framework of rainbow gravity, we need to have w = —^ or f{po) —>■ oo corresponding to po = (|^) 
The perturbations in the cosmic scale factor a{t) and the energy density p(t) can be written as 


Substituting (1551) into the equation 
equation 


a{t) ao(l + Sa{t)), 

P{t) ^ Poi^ + Spit)). (33) 

after linearizing the perturbation terms, we obtain the following 


2k ^ SttG , 

2 OG — PO^P ’ 

% 3 


(34) 


This indicates, with respect to the positiveness of k, that the sign of variation of the scale factor must 
be opposite to the sign of variation of the matter density. We must also apply the same method on 
the equation (1241) . however before doing this let us introduce the following replacement due to the 
perturbation in the rainbow function f{p) 


r\p) 


f {Po)--^^P$Sp 


such that 


r^ipo) = 1 - -epj 


(35) 


(36) 


Also, by combining the energy conservation equation (I25|) with the first rainbow modihed FRW equation 
we obtain 


d f ddf p 3(l + 3a;) if 
-: =-^PM 


a f a dp' f 4 

Then, after applying perturbation (l33l) . we have 


SttG 


a I 

af 


3(1 + 3w) 


1 




^Po (1 + / (Po) - nPo^P 


{ ^^Po(l + M “ 2(5a) 

V 3 


(37) 


(38) 


By substituting (1511) and (1551) into the above expression, one finds that --4=0. It seems that this 
result is independent of any specific form of rainbow function /. Therefore, it can be seen that the last 
term of the second rainbow modified FRW equation (l24l) does not contribute in derivation of stability 
conditions of Einstein static Universe in the framework of DGR scenario. Now, putting perturbation 
equations (l33l) into (l24l) and using (1551) . we get 


Sd = — 


47rG, 


(1 + 3a;) /o ^(po) - -^Po Po^P . 


(39) 


where inserting the equation (1541) into (1591) and neglecting the non-linear perturbation terms, results in 
the following differential equation 


10 1 

Sd -2 (1 + 3a;) ( 1- tt^Po ) = 0. 


It is obvious that, for the cases ^ 0, i.e. epi —>• oo. Eg. (1501) will take the following form 

k 

Sd - 2 (1 + 3oj)Sa = 0 , 

Og 


(40) 


(41) 


which refers to the oscillatory modes of ESU in the framework of standard GR for a; < —1/3. In order 
to have the general oscillating perturbation modes in the framework of the DGR scenario, the following 
condition should be satished 

- + >0, (42) 
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which results in the following solution for the equation (HOI) 


5a = Qfie*'*'”* + 026 


(43) 


where oi and 02 are integration constants and 70 refers to the frequency of oscillation around the stable 
ESU as 

Now, we can analysis the stability condition (1421) . Eas. (l3T]) and (1321) will play the role of two fundamental 
constraints to achieve our goal. There are two possibility in order to satisfy the equation (1321) : the first 
one is that a; = — ^ and the second one is that the matter density of the ESU be po = (|^) . If we 

set w = — i, then the condition (I42|) is automatically violated and we obtain (5d = 0 implying that there 
is no stable ESU against the linear scalar perturbations in the form of relations (1551) . As it is clear 
from Eg. lHTl) . in the framework of the standard GR for w = — there is no stable ESU. For the second 

possibility as po = (|0 with respect to the positivity of the spatial curvature, due to the positive 
energy condition through Equation (l3T]l . we need that the barotropic EOS parameter satisfies w > —1/3. 
Overall, one can find that there is a stable ESU versus homogenous and linear scalar perturbations in 
context of the gravity’s rainbow scenario with choosing the rainbow function (1271) . In this case, we need 
two bounds which the first one in on EOS parameter as w > —1/3 and the second one is on the energy 
density of ESU as po < There are no such similar results in the framework of the standard GR 
theory. Unlike GR, here there is a stable solution against small scalar perturbations for a closed universe 
filled by the matter fields respecting the energy conditions. Also, within standard GR theory and even 
in many of the modified gravity theories, there is no such upper bounds on energy density of ESU po- 
In the following, we want to investigate the stability conditions of Einstein static Universe against the 
linear homogeneous scalar perturbations in the context of DGR, in terms of another rainbow function 
proposed by Ling and et al in [35], which is evolving with cosmic time as 


/=/ru57. 


(45) 


Gonsidering previous arguments after placing the average energy e = |cp4 m, it can be rewritten in 
terms of energy density p as 


f{p) 



(46) 


where x = is a constant. Let us point out that rainbow function (1461) has application in the 
framework of black hole physics leading to valuable phenomenological outcomes, e.g. see |3SJ [5S] . 
Moreover, it also as previous rainbow function has a theoretical trait. By reaching the energy density p 
or e to its largest value, /(p) in (l46l) becomes zero which leads to vanishing spatial-like components of the 
rainbow metric. This means that FRW metric modified by (l46|) or (|45|) . is a smooth and differentiable 
metric but is not invertible. Moreover, by applying the rainbow functions (|45l) or (I46|) . we will be faced 
with a degenerate metric at energy levels close to the Planck energy scale. Also, based on the fact 
that there is no observational evidence for lightlike dimension in nature, there might be a suppression 
mechanism for appearance of lightlike dimension at the quantum level of Universe [32] . 

As before, we first examine the Big-Bang singularity problem for such a choice of rainbow function 
within a flat rainbow FRW model. In the same way, for the rainbow function (I46L we get 


p = -4 



(47) 


which has only one fixed point p/^ = 0 and is not bounded, see Figure 2. As it is seen from the figure, 
the function 0{p) begins from the fixed point pf^ and ends with a singularity. So we must separately 
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calculate the time needed from an infinite to a finite value p* as well as the time needed from p* to the 
fixed point =0, i.e, 


and 



(48) 


(49) 


Here, we are dealing with a different situation compared to the previous one. We found that the integral 
(H51l converges which provide p* > 2 ^^ while the integral does not converge anyway on the given 
interval. This means that from p* to fixed point pf^ = 0 an infinite time is needed, i.e |t| 00 , while 

to get an infinite to a finite value p* it takes a finite time. Because the rainbow FRW metric includes a 
natural cutoff as the energy Planck, the constraint obtained for integral (H51) should lie in the interval 
< P* < Then, the solution derived from the rainbow function (I46|) for the flat rainbow 
Universe model can not be free of finite time singularity. We conclude that even assuming that integral 



Figure 2: The behavior of 0(p) (I47II in terms of p for the flat rainbow FRW Universe model. To simplify, we 
set SttG = 3 and y = 1/2. 


(l4^ results in an infinite time, the Big-Bang singularity issue is not canceled, yet. Because according 
to the discussion given in the previous case, one can realize that p/^ is a future fixed point and is not a 
past one. 

In the following, we study the DGR theory from the emergent universe point of view and investigate 
the required conditions for the stability of an ESU with respect to scalar perturbations, in the presence 
of rainbow function (1461) . By applying the scalar perturbation terms on the rainbow function (1461) . it 
takes the form 

^ f~^iPo)+ ^XPo^P , (50) 

where 

/"^(Po) = 1 + yXPl • (51) 

Now, inserting the equations (|50)) . (IbTl) and the perturbation equations (|3^ into the second rainbow 
modified Friedmann equation (1241) . we have 

Sa = --^(1 + 3w) ^(po + l)/"^(po) + PoSp , (52) 

where similar to the previous analysis, we neglected non-linear terms. Finally, by applying the equations 
(15^ and (IMl) . the above differential equation takes the following form 

6a-^{l + 3uj) (^1 + ^XPo'^ ^a = 0 . (53) 
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As can be seen, for the limit % —0, the above equation reduces to the oscillatory modes (HTl) . Then, 
in order to have a stable ESU in the framework of DGR scenario, the following condition should be 
satisfied 

- ^(1+ 3w) (l + > 0 , (54) 

where the frequency of oscillatory modes 70 for the rainbow function (|46ll reads as 

7o =-^(l + 3w) + ■ (55) 

It is seen that for ui = — i, the condition (IMl) is automatically violated and we get (5a = 0 through 
the equation (15^ . which indicates that there is no stable ESU against the linear scalar perturbations 
for this rainbow function. For the second possibility as po = (^x) with respect to the positivity 
of the spatial curvature due to the positive energy condition through equation (I31|) . we need that the 
barotropic EOS parameter w satisfies the condition w < —1/3. Therefore, here we have a stable ESU for 
the baratropic EOS parameter w < —1/3, denoting the phantom matter fields, with the bounded energy 
density of ESU po as po < As mentioned before, the existence of such constraint on energy density 
of ESU may be absent in other modified gravity theories. Here, such a bound on initial density po is 
arising from the energy dependent metric in gravity’s rainbow proposal which is one of the common cut 
offs of quantum gravity theories. Finally, it should be mentioned that choosing the rainbow function 
(l46t leads to the same result with GR so that a closed universe hlled by usual non-relativistic matter 
helds, is not stable against small linear scalar perturbations. 


3.2 In the Presence of Cosmological Constant A 

In this subsection, we want to consider the possible modification by the cosmological constant A(e) and 
investigate its effects on stability condition of an ESU. In presence of cosmological constant A(e), the 
second modified FRW equation remain without change. But, the first modified FRW equation (P51) 
and the conservation law (1251) . take the form of 


,a ^2 , ^ 1 SttGp 1 ^ A(e) 


(56) 


and 

p + 3p^(l + a;) = -Ag. (57) 

One may consider another rainbow function as h(e) so that A(e) = h(e)^A where A is the usual cosmo¬ 
logical constant. By this consideration, the modified FRW equation (l56l) will be 


, 0.2 k \ SnGp I h{e)^A 
^ ^72 ^ 3 72 3/2 


(58) 


By setting this new rainbow function as /i(e)^ = (1 -I- Ap), see second Ref. of [57] i for tfi® ESU described 
by (i = d = 0, we get the following equation from equation (1551) as 


k SttGpq , a 


(59) 


where A is a dimensional parameter. By keeping positive energy condition, the above equation implies 
that positive spatial curvature fc > 0, for two cases A > 0 and A < 0, is guaranteed only under the 
following constraints, respectively 

A > -— - — , (60) 

Pa Po 
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and 


( 61 ) 



IpaI Po ’ 


so that PA = and |/cia| = With note to the our previous result on upper bounds on energy 

density of ESU as po < the above constraints reads as follows 



(62) 


and 



(63) 


In DGR formalism of gravity, it is expected that parameter A to have the order of magnitude |A| ^ 

As it seems, this value of A be satisfied both the above constraints. Also, in the presence of cosmological 


constant, we recover Eq. (15^ from the second modified FRW equation (IMl) . Similar to the previous 
analysis, by inserting equation (1331) into the equation (155)) and neglecting the non-linear perturbation 
terms, we get 



(64) 


It is seen that for the case of positive spatial curvature fc > 0, the sing of variation of the scale factor 
5a is in contrast to the variation of matter density 5p, provided that 



(65) 


and 



( 66 ) 


where are related to the cases A > 0 and A < 0, respectively. For the case of positive cosmological 


constant, A > 0, by comparing the constraints (1601) and (1551) . one realizes that constraint (1551) can be 
satisfied via constraint (|65l) while its inverse is not true and so relation (I65p will be a tighter constraint. 
With the same reason, one can say that for the case of negative cosmological constant, A < 0, constraint 
(l6T|) is more tighter than (1551) . Finally, by using equation (1541) and considering equation (1331) into the 
second modified FRW equation (I24L after neglecting the non-linear perturbation terms, we arrive to 
following differential equation 



(67) 


By putting po = (|^) the above differential equation takes the form of 



( 68 ) 


In order to have a stable ESU against homogeneous and linear scalar perturbation described by the 
oscillatory modes, we have two possibilities as 


• For the case of A > 0 



I 



Pa’ 


or o,<-i 



(69) 


• For the case of A < 0 
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a; > —- 


A < 


IpaI’ 


or uj < — 


3’ 


A > 


IpaI’ 


(70) 


Now, in the same way as above by using the rainbow function (1271) . by choosing the rainbow function 
(ITO . we obtain the following differential equation 


ag V 1 + ^PA J 


1 + 3XP0 


(1 + 3uj)Sa = 0 , 


(71) 


where by inserting pg = (^x) ^ 1 coming from constraint (1321) for the rainbow function (I46L it can be 
rewritten as 


(5d — 


5k 


1 


(1 + 3uj)Sa = 0 . 


1 + ApA ^ 

This differential equation possesses the stable oscillatory modes under the following possibilities 


(72) 


• The case of A > 0 


1 ^ 1 

^ > -n’ A <-, 

3 PA 


1 , 1 

or UJ < A > -. 

3 PA 


(73) 


• The case of A < 0 


UJ > -- 



1 

or w < 



(74) 


We note that since the cosmological constant does not appear in constraint equation (15^ . the energy 
density of ESU is bounded as po < e^i or po < for the both rainbow functions (071) and (H51) . It is 
found that in order to have stable solutions for an ESU in the presence of a cosmological constant, beside 
the cutoff on energy density po, the conditions (|6^ . (ITO)) and (l73|), (174)) for each of rainbow functions (l27l) 
and (|46l) . should also be satisfied respectively. For the case of positive cosmological constant A > 0, by 
comparing conditions mentioned in (I69|) and ()73l) with the constraint ()65l) . one realize that with choice 
of rainbow functions ((771) and (H51) . the solutions with uj > 

are allowed, respectively. For the case of negative cosmological constant A < 0, constraint (EO , is more 
tighter than (1551) . In fact, we recall that the constraint (1511) is important in the sense that with respect 
to physical energy condition in quantum gravity regimes, ensures that fc > 0. So, it is clear that all 
the constraints on the parameter A in dZO]) and dUl) violate the constraint (ED. Overall, this means 
that in the presence of a quantum gravity modifications such as what is done by gravity’s rainbow, 
with respect to positive energy condition and positive spatial curvature fc > 0, which is one of the basic 
assumptions in modern version of emergent universe scenario, for the cosmological models possessing a 
negative cosmological constant, there is no stability for an ESU. 


4 Einstein Static Universe, Vector and Tensor Perturbations 

In the cosmological context, the vector perturbations of a perfect fluid having energy density p and 
barotropic pressure p = ujp are governed by the co-moving dimensionless vorticity defined as zua = aw. 
The vorticity modes satisfy the following propagation equation m 

w, + (l- 3cl)Hw^ = 0, (75) 

where = dp/dp and H are the sound speed and the Hubble parameter, respectively. This equation is 
valid in our treatment of Einstein static Universe in the framework of the rainbow gravity through the 
modified Friedmann equations (1751) and (1751) . For the Einstein static Universe with H = 0, the equation 
dZHD reduces to 

Wk = 0, (76) 
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where indicates that the initial vector perturbations remain frozen and consequently we have neutral 
stability against the vector perturbations. Tensor perturbations, namely gravitational-wave perturba¬ 
tions, of a perfect fluid is described by the co-moving dimensionless transverse-traceless shear tensor 
Sab = o-o’ab, whose modes satisfy the following equation 


Sk -I- Si/S^ + 



2k 


(1 -I- 3w)p — 2A(e) 
3 


Sk — 0, 


(77) 


where K. is the co-moving index —K^jd} in which is the covariant spatial Laplacian) [37] . For 

the Einstein static Universe {H = 0), this equation reduces to 


Kd 2k (1-I-3w)po ~ 2A(1-I-Apo) 


S. = 0. 


(78) 


Note that in the general tensor perturbation equation (EZD, the Hubble parameter H contains the 
rainbow function / through equation (1231) . But, because for the Einstein static universe we have H = 0, 
then the rainbow function / did not appear in the equation (1781) . Then, in order to have stable modes 
against the tensor perturbations, the following inequality should be satisfied 


K.^ ^ 2k ^ {1 + 3a;)po - 2A(1 -f Apo) 
% ag 3 


(79) 


Where, for a closed universe with fc = 1, considering the eigenvalue spectra K.^ 
n = 1, 2, 3,... [38], will takes the form of 

fd 2 ti 2 (1 T 3u^)po — 2A(1 -t- Apo) 

9 ^ n ■ 


n{n + 2) with 


(80) 


This inequality gives a restriction on the scale factor of an ESU in terms of its matter density and 
background cosmological constant. Even though the existence of open and flat solutions are possible 
in some modified gravity models, in the present modified gravity model it is forbidden by regarding 
the weak energy condition. Therefore, the stability analysis in the present model is restricted to the 
physically viable closed cosmological model. For the the eigenvalue spectra Kd for open and flat models, 
one is referred to [55] . 


5 Summery and Conclusion 

According to a modern version of emergent cosmological scenario proposed by Ellis and Maartens 
m, the early Universe before passing to the inflationary phase, has experienced an eternal Einstein 
static state rather than a Big-Bang singularity. More exactly, it refers to a static closed space in 
the asymptotic past (before entering the Universe to the period of inflation) known as “Einstein static 
Universe” (ESU). The initial conditions as the quantum gravity effects in the early high energy Universe, 
will influence the stability of this static state . For this reason, in the present work, we have examined 
the effects of an effective approach to quantum gravity proposed by Magueijo and Smolin [10], known 
as ’’gravity’s rainbow” or ’’doubly general relativity (DGR)”, on the stability of the Einstein static 
state against the linear homogeneous scalar, vector and tensor perturbations. In order to following our 
aim in the framework of DGR, we have needed to introduce appropriate rainbow functions. Then, we 
have considered two appropriate well known rainbow functions (1271) and (1461) . First, by following the 
phase space mechanism represented in |31j . we have investigated the fixed points belonging to a flat 
FRW model which is modified by the presence of the rainbow function (l27l) . Indeed, these fixed points 
are de Sitter space solutions of typical flat rainbow FRW model. Our results indicate that although 
the needed time to get a fixed point is infinite, it may not lead to the elimination of the initial finite 
time singularities (Big-Bang). This is because of the fact that the fixed points corresponding to past 
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singularities (Big-Bang) are unstable and may be collapse. Then, we studied the DGR theory in the 
context of the emergent universe scenario and tried to find its stable Einstein static universe with the 
required conditions. In order to find a stable ESU, as the first result, we realized that similar to standard 
GR scenario and independent of the choice of rainbow function, the positive spatial curvature A: > 0, is 
the only option respecting to positive energy condition po > 0 for an ESU in the framework of the DGR 
theory. Eor the case of zero cosmological constant A = 0, we found that in order to achieve a stable 
ESU against homogeneous and linear scalar perturbations, we need a conditions on the energy density 
of ESU and the barotropic equation of state as po = (|^) and w > — j, respectively. Similarly, for the 
rainbow function (H51) . we also verified that Big-Bang singularity, can still exists. Unlike the previous 
case, the time needed to get the fixed point is not infinite in this case. Furthermore, the only fixed 
point which is revealed in the presence of rainbow function (1461) . is a future fixed point and is not a 
past fixed point. Then, by looking at the DGR theory from the emergent universe point of view, we 
studied the stability of the Einstein static universe and its required conditions. For this case, a stable 
solution for ESU is guaranteed under these conditions po = (^x) and a; < — ^ describing the exotic 
matter fields. As it is seen, stability conditions po = (|^) ^ and po = (^x) ^ are direct result from the 
idea of an energy dependent metric in the framework of gravity’s rainbow scenario. These results are 
equivalent to an explicit cutoff on the energy density of an ESU as po < and po < respectively. 
Given that, ESU point out a initial static state (or static closed space) of universe before getting into 
inflationary phase, so the existence of such explicit cutoff (or upper bounds) on po could be interpreted 
as a result of initial dominate quantum gravity effects such as “gravity’s rainbow”. In the following of 
our analysis, we take an energy dependent cosmological constant via the introduction of new rainbow 
function as h{e)'^ = (1 -I- Xp). It is seen that the positive spatial curvature, fc > 0, dictates an opposite 
sing for da relative to energy density Sp through the equation ()64l) . We find that in order to have 
stable solutions for an ESU,beside the cutoff on energy density po, the conditions (l69ll . dTHll and (IT^ . 
(l7^ for each of rainbow functions (l27ll and (l46ll . should also be satisfied respectively. In particular, for 
the case of positive cosmological constant A > 0, by comparing conditions mentioned in (l69l) and (l7^ 
with the constraint (1551) . one realize that with choice of rainbow functions (l?7l) and (H51) . the solutions 
with w>—i, ^ allowed, respectively. For the case of negative 

cosmological constant A < 0, constraint (1511) . is more tighter than (1551) . In fact, we recall that the 
constraint (1611) is important in the sense that with respect to physical energy condition in quantum 
gravity regimes, ensures that A: > 0. So, it is clear that all the constraints on the parameter A in (l70l) 
and dH violate the constraint (EH). Overall, this means that in the presence of a quantum gravity 
modifications such as what is done by gravity’s rainbow, with respect to positive energy condition and 
positive spatial curvature k > 0, which is one of the basic assumptions in modern version of emergent 
universe scenario, for cosmological models with a negative cosmological constant, there is no stability 
for an ESU. Finally, we investigated the the stability of an ESU in the framework of DGR versus vector 
and tensor perturbations. It is found that there is a neutral stability against the vector perturbations. 
In order to have the stability against the tensor perturbations, the scale factor of an ESU is restricted 
by its matter density and background cosmological constant. 
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